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1 Introduction 

At a summer school held in fall 1991 at the PTB in Braunschweig, Bernhard 
Kramer proposed a book project that should encompass novel solid state research 
topics ranging from quantum transport to quantum chaos. This very book finally 
appeared at the end of 1997 under the title "Quantum transport and dissipa- 
tion" B. Compiled between the book covers is a series of topics, some of which 
are usually not connected with each other in present day's research. The main 
subject in the first chapter is coherent transport in disordered systems. The last 
chapter, on the other hand, dwells on concepts such as phase space, Wigner and 
Husimi functions, and alike. In any case, in the description of disordered systems 
we do not find many works that explicitly make use of phase space concepts. This 
is so, although the existence of a mapping between the Anderson model and the 
kicked rotor indeed suggests that methods employed in quantum chaos may 
advantageously be used to elucidate the physics at work in disordered quantum 
systems. 

There exists, however, another, almost obvious physical motivation to uti- 
lize the powerful phase space concepts to study disordered quantum systems: 
as a function of the disorder strength, the nature of the eigenstates changes 
from ballistic to localized behavior, containing possibly in-between a diffusive 
regime ||. In the case of ballistic transport it is appropriate to think in terms 
of plane waves which are occasionally scattered by the weak disorder potential. 
Then a momentum space description clearly imposes itself. For strong disorder, 
however, the natural physical space is real space. A unified description, being 
valid at arbitrary disorder strength, can thus be achieved by focusing on suitable 
(quantum) phase space concepts. 

2 Phase space methods in quantum mechanics 

When an attempt is made to represent a quantum state in phase space, there 
are two observations which deserve our attention. First, phase space combines 
conjugate quantum observables which brings into play a Heisenberg uncertainty 
relation. While in classical mechanics, the state of a point particle can be de- 
scribed in phase space with infinite precision, this no longer holds true in quan- 
tum mechanics. Therefore, we expect some difficulties if for quantum systems 
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one attempts to obtain perfect resolution in phase space. Second, the full infor- 
mation about a quantum state is contained already in the corresponding wave 
function which can either be expressed in position or momentum space represen- 
tation. Thus, the phase space representation of a quantum state cannot provide 
more information than the wave function itself. This information, however, may 
be presented in a more advantageous form, as we shall demonstrate next. 



2.1 The Wigner function 



As a first possibility to define a phase space representation of a quantum state 
possessing the wave function ip(x) we introduce the Wigner function 



W(x,k) = I dye iky i>*{x + y-)ip{x- 



(1) 



Here, and in the following we use the wave number k instead of the momentum 
p = hk. The definition (Q) is useful and appropriate as can be seen by considering 
the moments of position and momentum observables. Upon evaluation of the 
momentum integral, one readily finds that 



x n W(x,k) = dxx n \ip{x)\ 2 . 



dxdk 

With a little more effort one further can establish that 
dxdk , 



(2) 



2tt 



-k n W(x,k) = dx dy 



dk 
2^ 



d™ ■ 



ip*(x)ip(x - y) 



dxip*(x) ( —— ) ?p(x) . 



(3) 



A generalization to moments containing both position and momentum operators 
is possible if a Weyl ordering for the operators is respected (for further details 
see [|). 

To start, it is instructive to consider a few special quantum states of simple 
structure which also play a key role in the discussion of disordered systems. 
Let us begin with a plane wave, i.e. ip{x) = exp(ifcoa;)/v / 27r. Inserting this wave 
function into (Flj) one finds 



W(x,k) 



^ e i(fe-feo)i/ = fifr _ fco ) 



(4) 



The intermediate result emphasizes the fact that the off-diagonal contributions in 
(0), which are parameterized by the coordinate y, contain the information about 
the momentum. For a localized state, ip( x ) — ~ x o)i the Wigner function 
emerges as 

W(x, k) = e 2ik{x - Xa] S(x - xq) = 6{x - x ) . (5) 
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Next we consider a quantum state localized at two positions, i.e. ip( x ) = [$( x + 
a) + S(x — a)]/V2. Proceeding as before, one obtains the Wigner function 

W(x, k) = - [8(x + a) + S(x - a)} + cos(2ka)S(x) . (6) 

The first two terms describe the localization of the particle at x = ±a. In addi- 
tion, there occurs an oscillatory term in the middle between the two localization 
centers which accounts for the coherent superposition of two localized states. 
This simple example demonstrates that the Wigner function generally is not 
positive even though it could be treated as a phase space density in (0) and 
(p|). This characteristic feature is a direct consequence of the fact that by virtue 
of the definition (|l|) one attempts to define a phase space representation which 
allows for perfect localization in position and momentum as indicated by the re- 
sults (S) and (||) . This formal attempt to circumvent the Heisenberg uncertainty 
relation is generally paid for by negative parts of the Wigner function. 

Thus far, we have considered a particle on a continuous and infinitely ex- 
tended one-dimensional state space. The situation changes when we try to apply 
phase space concepts to the Anderson model of disordered systems. For such 
a lattice model the integral in ([]]) has to be replaced by a sum. Moreover, in 
numerical calculations, the system size has to be taken finite. This results in a 
phase space that is twisted to a torus so that additionally periodic boundary 
conditions appear in momentum space apart from those usually imposed in real 
space. Therefore, interference terms in momentum space analogous to those ap- 
pearing in real space as in (^|) occur also across the boundaries of the Brillouin 
zone. Such artifacts can also be understood as arising from the discretization of 
the integral in ([!]) due to the lattice structure of the model. In particular for 
the states located at the boundaries of the Brillouin zone, the wave number k is 
so large that a discretized version of the Fourier integral no longer represents a 
good approximation. The usage of Wigner functions for lattice models therefore 
becomes problematic. 



In Fig. 1(a) we present an example for the Wigner function of an eigenstate 
of the Anderson model at relatively low disorder. While a certain spatial local- 
ization has already set in, two dominant wave numbers ±fc can still be clearly 
recognized. At k = 0, interference effects analogous to those found in @ are 
visible while the two low ridges at larger wave numbers represent the artifacts 
due to the lattice structure of the Anderson model discussed in the previous 
paragraph. 



2.2 The Husimi function 

The above discussed shortcomings of the Wigner function can be circumvented 
by use of the so-called Husimi function. The latter one is obtained from the 
Wigner function by a Gaussian smearing according to 

p(x,k) — — fdx'dk' exp 



x>) 



l\2 



2a 2 



-2a 2 (k-k'f 



W(x',k') 



(7) 
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Fig. 1. (a) Wigner function and (b) Husimi function for an eigenstate of the Anderson 
model 



Even though one might assume that one discards information by transforming 
the Wigner function W(x,k) into the Husimi function p(x,k), this is actually 
not the case as one can show by inverting (Q) : The Wigner function can then be 
recast as 



W(x,k)= (8) 

dudv dx'dk' - ■ , , , 

-~ + i(x — x)u + i(fe — k)v 



2tt 2tt 



■ exp 



p(x',k') 



The existence of the integrals is guaranteed by the inherent asymptotic decay 
of the Husimi function provided that the integrals over x' and k' are evaluated 
first. 

From a physical point of view, the most important aspect of the Gaussian 
smearing introduced in (Q) is that the Heisenberg uncertainty relation is now 
accounted for in a natural way. The Husimi function will never localize a quantum 
state in phase space beyond the limits set by the requirement AxAk > 1/2. In 
return, the Husimi function yields a nonnegative- valued phase space distribution. 
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Indeed, by inserting (|l|) into (fjj) one can express the Husimi function in terms 
of the wave function tp(x) as the squared quantity 



p(x, k) 



dx'ijj(x') 



(2 7 ra 2 ) 1 /4 



exp 



(x' - xf 



— ikx' 



(9) 



which manifestly is nonnegative. According to (|{j), the Husimi function may be 
viewed as the projection of the wave function onto a minimal uncertainty state 
or, using the language of quantum optics, onto a coherent state. At this point, 
the width a is still a parameter which can be chosen freely. We will fix its value 
later on when applying these phase space concepts to the Anderson model. 

Yet another aspect of the Gaussian smearing pertains to the artifacts men- 
tioned in the previous section for the Wigner function on a lattice. While in (Q) 
all values of y contribute, this is no longer the case for the Husimi function where 
the convolution with a Gaussian restricts the difference in the arguments of the 
two wave functions. In fact, even the interference term appearing in the Wigner 
function (||) for a continuum model is strongly suppressed in the corresponding 
Husimi function provided the distance a is sufficiently larger than the width cr; 
this is evident from 



p(x,k) = 



1 



exp 



(x - a)' 



2a 2 



exp 



2 exp 



(a; + af 

,2 



.JL_ A_)cos(2fca) 



(10) 



A com parison of the Wigner function in Fig. 1(a) and the Husimi function in 
Fig. 1(b) demonstrates the effect of the Gaussian smearing. The spatial localiza- 
tion and the two mainly contributing wave numbers are the dominant features of 
the Husimi function. The interferences leading to negative contributions to the 
Wigner function have almost disappeared except for some wiggles which, how- 
ever, do not render the Husimi function negative. Finally, the features stemming 
from the lattice have now disappeared. 



2.3 Inverse participation ratio 

The Husimi function, in particular for quantum states at intermediate disorder 
strengths, does contain a rich structure. In general, it is desirable for various 
reasons to reduce this wealth of information and to characterize the state by a 
single number which is derived from the Husimi function. First, in clear contrast 



to the one-dimensional case depicted in Fig. 1(b) it is generally not possible to 
visualize the Husimi function for systems in two and higher space dimensions. 
Second, the details of the Husimi functions depend sensitively on the chosen 
disorder realization. In order to perform averages, it is necessary to reduce the 
characterization of the states to a number. 

Since the Husimi function is nonnegative, in principle all quantities familiar 
from classical phase space analysis can be employed in the quantum case as well. 
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One possibility is given by the so-called Wehrl entropy Q 

dxdk 



2tt 



p(x,k)ln\p(x,k)], (11) 



which has been used, e.g., in the discussion of the driven rotor g]. The above- 
mentioned existence of a mapping between the kicked rotor and the Anderson 
model has motivated a study of the latter by means of the Wehrl entropy . 

From a numerical point of view, it is advantageous to linearize the Wehrl 
entropy (|Tl|). Replacing ln(x) by its linear approximation x — 1, one obtains 
— xln(x) x — x 2 . The approximation shows qualitatively the some behavior as 
the original function and, in particular, yields the same values at the boundaries 
at x — and 1. Performing this linearization in (fill), the first order term gives 
rise to a constant due to normalization and we are left with the second order 
term, the inverse participation ratio (IPR) in phase space 

p = J^WM 2 - (12) 

While both the Wehrl entropy S and the phase space IPR P are entirely deter- 
mined in terms of the wave function ip(x), it is necessary for the evaluation of 
the Wehrl entropy to first calculate the Husimi function explicitly. This requires 
a significant numerical effort which can be reduced by resorting to the phase 
space IPR ||. In turn, this makes possible the study of the Anderson model 
in two and even in three dimensions ||. Inserting (||) into ( |l2|) one obtains the 
phase space IPR expressed directly in terms of the wave function as 



p = -jr7=- 1(1,1 



dv ip ( — - — I tp ( — — — I exp 



(13) 



Apart from the numerical aspects, another advantage relates to the fact that 
inverse participation ratios may be defined as well in real and momentum space, 
i.e., 

P x = fdx\^(x)\ 4 and P k = f dk |^(fc)| 4 , (14) 



where tp(k) is the wave function in momentum representation. The availability of 
the inverse participation ratio in different spaces allows for instructive physical 
comparisons. The real space IPR P x is a well studied quantity as it is related 
to the return probability of a diffusing particle (k} while the phase space IPR 
has recently been used to describe the complexity of quantum states [ fu|]l2]| . We 
add that a quantity containing information similar to the phase space IPR can 
be defined on the basis of marginal distributions in real and momentum space 
P). 

In order to demonstrate that it is justified to refer to the above-mentioned 
quantities as inverse participation ratios, we consider a few special cases for 
quantum states on a one-dimensional lattice with L lattice sites. Let us first 
assume that the state is localized on a single site. Then we find P x = 1 and 
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Pk = l/L. The inverse of these results indeed indicates that in real space the 
state occupies only one site while in momentum space L sites are occupied. In 
phase space, one obtains P = 1/2 y/ncr which primarily reflects the Gaussian 
width in real space. 

As a second example, we consider a real-valued ballistic wave function at 
momentum k, i.e. ip(x) = \/2jL cos(fcc) . In real and momentum space, the 
inverse participation ratios become P x = 3/2L and Pk = 1/2, respectively. The 
first result accounts for the nonuniform extension in real space while the latter 
indicates the contribution of the two momenta ±fc. For k ^ 0, ±7r, the phase 
space IPR P = (ay/n/L)[l + 2exp(— 4a 2 k 2 )] depends on the distance between 
the two contributing momenta. 

So far, we have not specified the width a entering both, the definition (^) 
of the Husimi function and the results for the phase space IPR P mentioned 
in the previous two paragraphs. Although one is rather free in choosing a, the 
most impartial choice consists in selecting an equal relative resolution in real 
and momentum space, i.e. Ax/L = Akjl-n. Together with the requirement of 
minimal uncertainty, AxAk — 1/2, one finds Ax = a = y/L/ 4tt. For our two 
examples mentioned above (which will appear later in the limits of vanishing 
and very strong disorder) this implies, that the IPR P will scale like 1/VL as 
a function of system size. This result for one-dimensional systems is found to 
generalize to P ~ L~ d / 2 in d dimensions. 

An intuition for the differences between inverse participation ratios in real 
and momentum space, on the one hand, and the IPR in phase space, on the 
other hand, can be obtained by considering the resolution provided by the wave 
function and the Husimi function in real and momentum space. The vertical 
gray stripe on the left in Fig. corresponds to the real space wave function. 
Independent of system size, this wave function always leads to perfect resolution 
in real space. The price to be paid, however, consists in the absence of any 
resolution in momentum space. The converse holds true for the momentum space 
wave function. The situation is quite different for the Husimi function symbolized 
by the gray disk in Fig. |[ Although the discussion can be generalized 
we will consider the case where a ~ \[L. Then, structures occurring on 
lattice points, or less, cannot be resolved. While this resolution becomes worse 
as the system size is increased, the resolution relative to the system size improves 
with 1/yL. For large system sizes, the phase space approach therefore allows for 
a detailed description of a quantum state both in real and momentum space, 
while a specific wave function representation will always neglect one of the two 
possibilities. 

3 Anderson model in phase space 

In the preceding section we have demonstrated that a phase space approach 
indeed presents advantages when both, real and momentum space properties are 
of interest as it is the case for disordered systems when the disorder strength 
is varied from zero to infinity. Although, in principle, the Husimi function is 
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Fig. 2. The different phase space resolution provided by real space wave function and 
Husimi function is visualized by the gray shaded areas on the left and right, respectively 



equivalent to the wave function, it will make the relevant information more 
readily accessible. We next want to illustrate this very fact by considering the 
Anderson model for disordered systems. Its Hamiltonian |l5j 

H = -t (|x')(x| + |x)(x'|) + ^5>„|x)(x| (15) 

<x,x'> n 

is defined on a d dimensional square lattice with L sites in each direction. In 
order to avoid boundary effects, periodic boundary conditions are imposed. The 
first term describes the kinetic energy which allows for hopping between near- 
est neighbor sites < x,x' >. In the following, the hopping matrix element will 
set the energy scale, i.e. t = 1. The second term on the right-hand side of ( |l5|) 
represents the disordered on-site potential where the energies v n are indepen- 
dently drawn from a box distribution on the interval [—1/2; 1/2]. The respective 
disorder strength is then determined by W. 



3.1 Husimi functions 



Before we investigate the inverse participation ratio for the Anderson model, it 
is instructive to first take a look at the underlying Husimi functions. Since a 
visualization is readily possible only for one-dimensional models, we compare in 
Fig. H the one-dimensional Anderson model and the Aubry- Andre model |Tq| . 
The latter is based on a periodic potential A J2 n cos(27r/3n)|n)(n| replacing the 
disorder potential in (|l5|). Choosing (3 as the golden mean, (V5 - l)/2, the 
potential is incommensurate with the underlying lattice and displays a phase 
transition from delocalized to localized states [|l6) . We have chosen this model for 
comparison since its phase space IPR is comparable to that obtained for the two- 
and three-dimensional Anderson model |jjfl . In Fig. g, the Aubry- Andre model 
therefore serves as a substitute for the higher-dimensional Anderson models. 

Before we discuss in further detail the Husimi functions depicted in Fig. [?], 
we mention that for real wave functions the Husimi function, according to its 
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Anderson model Aubry- Andre model 




Fig. 3. Husimi functions for (a— e) the Anderson model with disorder strength 
W = 0.1,1,2,4, and 40 and (f— j) the Aubry-Andre model with potential strength 
A = 0.1,1.5,2,2.5, and 10 for a state close to the band center and a system size 
L — 377. In view of the symmetry p(x, —k) = p(x, k) only the upper half of the phase 
space is shown 
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definition (|^), respects the symmetry p(x, k) — p(x, —k). It is therefore sufficient 
to plot the Husimi functions solely for positive momenta from which the lower 
part may be reconstructed as a mirror image with respect to k = 0. 

The Husimi functions for the limiting cases of very weak and very strong 
potential depicted in Figs. 1(a) and 1(e) for the Anderson model and Figs. ||(f ) 
and ||(j) for the Aubry- Andre model can readily be interpreted in terms of states 
localized either in momentum or in real space. In addition, in Fig. 1(f) the 
coupling between plane waves due to the periodic potential becomes visible while 
the localization at two sites seen in Fig. ||(j) arises because the calculation was 
restricted to the class of antisymmetric states fl^ . Furthermore, we keep the 
state number fixed while passing through avoided crossings so that states may 
be localized at different positions as a function of potential strength (cf. e.g. 
Figs. 1(d) and 1(e)). 

More interesting than the interpretation of the Husimi function in the limiting 
cases is a discussion of the transition from plane waves for a weak potential to 
spatially localized states for a strong potential. The one-dimensional Anderson 
model and the Aubry- Andre model display two different scenarios with the latter 
being also characteristic for higher-dimensional Anderson models (cf. Sect. 3.2). 

For the one-dimensional Anderson model, the Husimi function for a plane 
wave (Fig. 1(a)) contracts in real space (Fig. 1(b)) and displays a very well 
localized core at intermediate disorder strength W (Fig. 1(c)). A further increase 
of the potential strength leads to a spreading in momentum direction (Fig. 1(d)) 
and finally a state localized in real space (Fig. 1(e)) is approached. 

For the Aubry- Andre model the scenario is quite different. As mentioned 
above, already a weak potential strength A can lead to a coupling between plane 
waves of very different momenta (Fig. |3|(f ) ) . This type of coupling is accompanied 
by an increased filling of phase space (Fig. 1(g)) up to potential strengths A 
located just below the critical value of A = 2 where the localization transition 
takes place (Fig. 1(h)). There, for increasing system size, a more and more abrupt 
contraction in phase space occurs, cf. Fig.|3|(i). Finally, for very strong potentials, 
the Husimi function of a localized antisymmetric state depicted in Fig. ||(j) forms. 



3.2 Inverse participation ratios 

While the Husimi functions depicted in Fig. 1 nicely elucidate the phase space 
behavior of a state as a function of potential strength, a more quantitative anal- 
ysis is desirable. On the basis of inverse participation ratios, a study of the 
higher-dimensional Anderson model becomes feasible and averages over disorder 
realizations can be performed. This in turn will allow us to draw general con- 
clusions about the physics that rules the localization transition in the Anderson 
model. 

In Fig. |i| we present distributions of the inverse participation ratios P x , P, 
and Pk in real space, phase space, and momentum space, respectively, for the 
Anderson model in one, two, and three dimensions. The distributions have been 
obtained by diagonalizing the Anderson model for 50 disorder realizations in 
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Fig. 4. Distribution of inverse participation ratios in real space (left), phase space 
(middle), and momentum space (right) for the Anderson model in d = 1 (upper row, 
L = 2048), d = 2 (middle row, L = 64), and d = 3 (lower row, L = 20) 



d = 1 and 2 and taking half of the eigenstates around the bandcenter for each 
realization. In d = 3 it was sufficient to consider only 20 disorder realizations. 

Upon comparing the real space IPRs for the three different dimensions, no 
qualitative differences can be seen. The overall behavior shows an increase of P x 
(left column) with increasing potential strength reflecting the spatial localization 
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of the wave functions. Correspondingly, the momentum space IPR (right column) 
decreases with increasing potential strength. 

In clear contrast, the phase space IPR behaves qualitatively different for the 
Anderson model in d — 1 and d > 2. For d = 1, the peak in P is consistent with 
our observations for the Husimi function in Figs. [|(a-e), namely the contraction 
of the Husimi function at intermediate disorder strength. 

For d — 2 and 3, the behavior of the phase space IPR rather corresponds 
to what we have observed for the Husimi function of the Aubry- Andre model 
in Figs. ||(f-j). With increasing potential strength, the Husimi function spreads 
in phase space implying a decrease of P. Then, at the transition, a sudden 
contraction occurs, implying a jump to relatively large values of P. This scenario 
can be understood in physical terms as follows. The large spread in phase space is 
associated with the presence of a diffusive regime which is known to exist in d = 2 
and larger. This interpretation is corroborated by a comparison with results 
from energy level statistics JM. The jump in P, on the other hand, indicates 
the Anderson transition Jl^JlJj]. Here, one might object that the case d — 2 is 
the marginal case where strictly speaking no phase transition occurs. One finds 
indeed, that both, the minimum and the maximum of P, shift towards vanishing 
disorder as the system size is increased [|o| . In contrast, in d = 3, the minimum of 
P shifts to larger disorder strength W , while the maximum shifts in the opposite 
direction, in agreement with the existence of the Anderson transition. 

A comparison of the Husimi functions for the one-dimensional Anderson 
model and the Aubry- Andre model reveals the differences between the Ander- 
son model in d = 1 and in higher dimensions. In the one-dimensional Anderson 
model, a weak disorder potential predominantly couples a plane wave to other, 
energetically almost degenerate plane waves. Therefore, the coupling occurs be- 
tween plane waves of almost the same momentum. Due to the finite resolution 
in phase space, this coupling cannot be observed in momentum space. In real 
space, however, the coupling results in large scale variations of the Husimi func- 
tion which in the end cause a contraction in phase space as discussed above. In 
clear contrast, for the Aubry- Andre model the situation is quite different: Here, 
the coupling predominantly occurs to distant momentum values as can be seen 
in Fig.j3|(f). This finally leads to a spreading of the Husimi function in phase 
space [P-4| . The very same scenario is valid for the Anderson model in two and 
higher dimensions. Again, there exist energetically almost degenerate momenta 
which are very different from those present in the original plane wave. The cou- 
pling due to a weak disorder potential thus again leads to a spreading in phase 
space. On the other hand, perturbation theory for strong disorder shows that 
the limiting value of the phase space IPR for large W is approached from above 
[ po| . Therefore, a jump from small to large values of P arises and thus a phase 
transition (in d = 3) is to be expected at an intermediate disorder strength W . 

In conclusion, these considerations in phase space in terms of concepts such 
as the Husimi function and the corresponding inverse participation ratio prove 
indeed very valuable in order to explore in greater detail the physics for a class 
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of quantum systems which is so dear to Bernhard Kramer, namely disordered 
quantum systems [||. 
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